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Abstract 


The transport of galactic cosmic-ray helium nuclei and their 
secondaries through bulk shielding is considered using the straight-ahead 
approximation to the Boltzmann equation. A data base for nuclear 
interaction cross sections and secondary particle energy spectra for high- 
energy light-ion breakup is presented. The importance of the light ions 
H, H. and ' He for cosmic-ray risk estimation is discussed, and the 
estimates of the fractional contribution to the neutron flux from helium 
interactions compared with other particle interactions are presented us- 
ing a 1977 solar minimum cosmic-ray spectrum. 


Introduction 


The description of the transport of solar and cosmic radiation in spacecraft and satellites is 
important for estimating radiation biological harm to astronauts. The galactic cosmic rays (GCR’s) 
contain a large alpha-particle component that consists of approximately 10 percent of the total flux; 
this alpha-particle component is second in size only to the hydrogen component (ref. 1). Light-ion 
transport will be important for studying the large solar particle eyents that contain a significant 
alpha-particle component in addition to the dominant protons (refs. 2 and 3). Target fragmentation 
effects also produce a large fluence of alpha particles and other light, ions; these particles and ions 
are expected to be biologically damaging because their spectra peak at low energies. An accurate 
physical description of the transport alpha particles is thus needed for estimating the radiation fields 
to be encountered by astronauts in spaceflight. The production of 2 H and 3 He from alpha-particle 
fragmentation is also useful to estimate the natural abundance of these elements in the GCR’s. 


For high energy, charged-particle transport, the straight-ahead approximation is accurate The 
solution of the Boltzmann equation for particle propagation in shielding materials has been developed 

no™™VT d r Lamkm ' Fef ' 4 l a " d Wilson ( rcf - 5 )> ancl a s P ace radiation transport code called 
HU Y IN RN (rels. 6 and 7) provides an accurate numerical procedure. The accuracy of the solution 
is also dependent on the nuclear interaction cross sections for fragmentation reactions in shielding 
materials (ref. 8). Presently, there is a scarcity of experimental measurements for fragmentation cross 
sections; therefore, theoretical models must be used to provide the large data base of cross sections for 
transport in spacecraft shielding and tissue. One approximation made for heavy projectile fragments 
is that the secondary spectrum is a delta function corresponding to the beam velocity that closely 
resembles experimental observation; however, for light secondaries, the spectrum is much broader. 
This broader spectrum places an extra burden on the development of a data base because the 
secondary spectrum, as well as the total yield of each species, must be represented in the code. 

In this paper, we discuss the development of a parametric data base for alpha-particle interactions 
and secondaries He ( h ), H (£), and “H (d) in arbitrary materials based on theoretical models and 
muted experimental information. A description of the extension of the BRYNTRN code for light-ion 
transport is also described, and predictions for cosmic-rav helium transport are discussed. Alpha 
particles will be the dominant source of high-energy neutrons (>100 MeV) in the GCR’s and an 
important source of low-energy neutrons. Estimates of the fractional contribution of secondary 
neutrons compared with neutron production from other particle types will also be discussed Light- 
ion secondaries from target fragmentation are included in the transport equations. In previous 
treatments (refs. 6 and 7), these secondaries were included as a local source. The new procedure has 
the advantage of including the effects of any further nuclear reactions and their contributions to the 
linear energy transfer spectrum. 



Transport Equations and Solutions 

The propagation of high-energy 4 He and its secondaries through bulk matter is described by the 
Boltzmann equation, which in the straight-ahead approximation is of the following form (refs. 5 
and 6): 


ic ~ V i§E S[E) + aj{E) } ^ E) = ^ if fjk (£ ’ E>) * k (x ’ E>) dE ' 


(1) 


where Uj (which denotes the range scaling parameter) is equal to Z^/Aj. The terms Z and A are the 
charge and the mass number, respectively. In equation (1), S{E) is the proton stopping power, trAE) 
is the total cross section, <pj{x, E) is the differential flux spectrum of type j ions, and f jk {E, E ) is 
a differential energy cross section for the redistribution of particle type and energy. 


Utilizing the definitions 


r E dE ' 

'o sW) 


( 2 ) 


* j (x > r) = S{E l )<t> j (x,E) 


(3) 


and 

f jk 0 r,r’)=S{E') f jk (E,E ') 
allows equation (1) to be rewritten as 


d d ^ 

— i ”77 — \ & i 

dx 3 dr 3 



Vj{x,r) 


^j r fjk(r,r)V k (xy) dr' 


(4) 

(5) 


This equation is solved by Wilson (ref. 5) and Wilson et al. (ref. 6) as 

<Iq (X, r) = exp(— Cj(r, x)\*j (0, r + v jX ) + ^ ^ £ exp[-Cj(r, x))f jk (r + v 3 z, r') (* - 

where the exponential quantity is the integrating factor given by 

(j(r,t) = ajir + Vjt') dt ' 


2 , r') dr dz 


( 6 ) 


(7) 


Equation (7) is the basis of the numerical procedures for propagation of the solution at ’fj ( x,r ) 
to ^ j (x + h,r). If we choose h to be small such that 

aj (r) h 1 

the solution at x + h (ref. 6) is approximately 

tyj(x + h,r) = exp[-Cj(r, h)} Vjixj, r + vjh) 

+ £ exP {- [<* ^ ( r ’ £)] } T ( r ’ "A Tk 


+ 


v j h 


x f k V k {x,r')dr' 

•' r k+(i'jh/2)+(i'i c h/2) 


(9) 
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for charged-particle propagation and 


^n{x + h,r) S exp [— Cn(r, h)] <S> n (x,r) 

+ ?f^H4’l) + *H)]} 

x fjk ( r > r ' ) (z, r' + v k dr ' (10) 

for neutron propagation. In equation (9), F is related to the cumulative spectrum F as 


Fjj (h,i 


f h - 

J {) fij ( r + 2 ’ r ') dz = Fij (r + h, r') - F v] (r, r') 


(ID 


with 


Fi 


re(r) 

i( r ’ r ) = j Q fij (E, E') dE 


( 12 ) 

where e(r) is the energy associated with the residual range r and E' = sir'). The numerical 
procedures for evaluating equations (6) to (12) have been described in references 6 and 9. 

For nuclear secondaries produced from target nuclei with A > 4, the probabilities for nuclear 
collisions are small because these will be low-energy ions (e.g., typically only a few McV). The 
differential flux for these ions can then be solved in closed form for the target fragments with A > 4 

cliS 

1 poo 

1 V / 


<pj(x, E) = 


j{E) ^J E 


fkj ( f \ E') <p k (x, E) dE' 


(13) 


The light-ion target fragments have been treated in the past (refs. 6 and 7) by using equation (13). 
however, here they will be included in the transport solutions. 

^5??TRM dCS ? rib p thC nUmCriCal re P rescntation of Particle ranges and stopping powers in 
the UKYNTRN code. For energies greater than a few MeV/amu, Bethe’s theory using the Born 
apprmamation is adequate if the appropriate corrections to Bragg’s rule, the shell corrections, and 
e effective charge are included. Proton stopping powers are taken from the parametric expressions 

0 ndersen and Ziegler (ref. 10). A modification to their shell corrections, however, has been added 
to ensure a smooth transition to Bethe’s asymptotic formula. For alpha particles, the electronic- 
stopping power cannot be derived from the proton stopping power at low energies because of the 
neglect of higher order Bom terms. Instead, we use Ziegler’s method (refs. 11 and 12) based on fits 

o the experimental data For low energies, the nuclear stopping theory used herein is a modified 
form of the theory of Lindhard et al. (ref. 13). 

Light Ion Fragmentation on Nuclear Targets 

Although the alpha particle is the most compact and tightly bound nucleus, this does not preclude 
its breakup in the held of target nuclei. Fragmentation of cosmic-ray alpha particles may lead to 
an important source of high-energy neutrons that will have large ranges and high values for relative 
biological effectiveness. Theoretical models of light-ion cross sections have been developed (refs. 14 
to zU), and we now consider these models for providing a data base. 

The fragmentation of 4 He is simpler than that of a heavier nucleus because only a small number 

01 final states can occur. These reactions are 

4 He + AT 
3 He + n + X 
4 H + p -T X . 

2 H+2 H + X } (14) 

2 H + n+p + AT 

n+n+p+p+X 


a + T 
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where T is the target and X is the final target. Each of the reactions in equation (14) can occur with 
or without meson production if sufficient energy is available. The reactions m equation (14) are not 
exhaustive of the absorption processes; most notable of these reactions are the compound nuclear 
and pickup channels that are important at low energies. A model for the tw<>body dissociation of 
light ions has been developed (refs. 16, 17, and 20), and this model describes the first three reactions 
in equation (14), which we now discuss. The fragmentation of the composite secondaries of He ( H, 
3 H, or :i He) is described by the same model. 


For an inclusive reaction involving the two-body dissociation of the projectile P, we write 


P +T -> a + b+ X 


(15) 


where a and b are assumed to be the clusters that are initially present in the projectile and X is t re 
final unobserved target state. We consider the case in which a is the observed projectile fragment in 
the measurement. Note that we must sum the unobserved target states to evaluate the cross section. 
We should also consider summing possible states of the particle b. Using relativistic kinematics, the 
transition matrix for equation (15) is related to the Lorentz invariant momentum distribution for 
producing the fragment a by 



where f and i label the final and initial states, 0 is the relative projectile-target velocity, and tj, 
is the transition matrix. The summation in equation (16) is over the possible configurations of the 
target in the final state. 

The transition matrix can be written as a three-body problem of aT, bT, and ab interactions 
when rearrangement channels are neglected and with the understanding that all target final and 
intermediate states must be summed. Using the Faddeev method, we consider the multiple-scattering 
series generated by the coupled set of integral equations 



T = r° +T b + T 1 

(17) 

with 

r a = T b r f + T b r (j h + ^ 7 ) 

(18) 


T b = f aT + T aT Go (?" + r T ) 

(19) 


T 1 = T ab + T ab Go (r" + T 6 ) 

(20) 


where T r Tvr, and T ab are the “two-body” amplitudes that are the transition operators for aT, 
bT, or ab scattering, respectively, in the projectile-target Hilbert space. Solutions for the amplitudes 
related to f aT and T br arc found (refs. 15 and 19) from Watson’s multiple-scattering series in 
the eikonal approximation. The amplitude related to T ab will only be important at small ab relative 
momentum, and it is evaluated in a separable potential model. The Green’s function G 0 is evaluated 
in the impulse approximation that neglects the binding of the nuclei of a, b, and T with respect to 
their kinetic energies. 

We consider the leading terms of equation (17) by truncating as 

T=(l + f ab Go) (faT + f b r + f aT G 0 f bT + T br Gi)T a ’j^ (21) 

which considers the scattering to all orders but neglects reflection terms between projectile and target 
nucleons. This equation assumes that the final-state interaction (FSI) between projectile fragments 
occurs only after interactions with the target. 
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The evaluation of the momentum distribution, equation (16), using equation (21), proceeds 
after an mtroducUon of the vertex functions for the projectile decay and the inclusive scattering 
distributions for a or b interacting on the target. The vertex function is parameterized as a sum of 
Yukawa functions with the proper asymptotic behavior at large ab separation distances. Reference 20 
gives the values for parameters of the Yukawa representation of the vertex functions for light ions. 

The inclusive scattering distributions for a or b reacting on the target are found (refs. 18 and 19) 


<i 2 < t 

dn dE 


K 2 

(2?r) ; 


/ At 

d b d 2 b f cxp['iq(b - b')] exp{i[*(b) - * f (b)]} ^ \V m (b ? b', 


( 22 ) 


m=l 


where K is the relative wave number, y is a mean distorting wave in the eikonal approximation 
q is the momentum transfer, u> is the energy loss of the projectile fragment, and W m are the 
collision terms for the mth order inelastic scattering. The collision terms are related to the two- 
body amplitudes, and the form factors of the projectile cluster (a or b) and the ground-state response 
functions (Greens functions) of the target (refs. 18 and 19). The elastic fragmentation terms arc 
included using a hrst-order optical model solution for the aT or bT amplitude. The on-shell spectra 
represented by equation (22) are used to approximate the off-shell amplitudes f aT in equation (21). 

igures 1 and 2 show the calculations of total a-nucleus reactions as a function of the invariant 
momentum transfer for several targets compared with the data of reference 21. Calculations of the 
inelastic part arc from equation (22) after integrating the energy loss. The elastic part is calculated 
in the hrst-order optical model. The results indicate the accuracy of using equation (22) for inclusive 
scattering at high energies. 

Figure 3 shows the angular distribution for 3 He production in a + *H reactions at 1 04 GcV 
compared with the data of reference 22. The calculations show an important interference effect 
among the terms m equation (21). Figure 4 shows the calculations of the invariant momentum 
distribution in the longitudinal direction near 2.0 GeV incident energy with data from reference 23 
and figure 5 gives data from references 24 and 25 for the transverse direction. Theoretical calculations 
provide a good representation of the data and indicate the importance of final-state interactions 
Results at other energies as described in reference 20 indicate a similar success for the light-ion 
breakup model just described. 


The spectral distributions of secondaries used for the solution of the Boltzmann equation (5) are 
obtained from equation (14) after integrating the solid angle. Total fragmentation cross sections 
are obtained from equation (14) after integrating both the solid angle and the momentum of the 
fragment. Table 1 presents the results of calculations for 3 H and 3 Hc production compared with 
experimental data (refs. 22, 26, and 27). 


Parameterizations of Interaction Cross Sections 

\\c next discuss parameterizations of interaction cross sections for light ions in common shielding 
materials. The work of Meyer (ref. 26) gives a complete summary of a- 1 H cross sections based 

qooTt^V 011101118 takCn Untl * 1972 ' Parameteri7 -ations of 3 H and 3 He production of 'll below 
300 MeV/amu are discussed in reference 28. An extensive list of earlier references of experiments is 
given in reference 26 Reference 29 discusses more recent experiments in absorption cross sections 
between 18 and 48 Me V. Reference 30 describes deutcron production at 1.4 GeV/amu which was 
measured with the inclusive deutcron production cross section reported at 30.64 ± 0.62 mb The 
results for the A = 3 fragments are given in references 22 and 27, and these results postdate the 
compilation by Meyer. The most important shortcomings of the data base for ‘‘He- 1 !! interactions are 

ngh-energy measurements above a few GeV and a complete absence of data for nucleon production 
cross sections. 
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Using our theoretical estimates and the existing data, we parameterize the fragmentation cross 
sections for 3 He, 3 H, and 2 H production of X H as 


^He = 425 


1 + exp [(E t h - E) / 6.8] 
0.51 


- 1 


as H = 15.5 


1 + 6.7exp(-.E/34) 

2 

[ 1 + exp [(£</, - E ) /7] 
0.45 13 


1 - 


1 + 7 exp (-E/ 55) 


E 


1 + 0.36 \l J exp [- (E - 780) /2300] (23) 

52U J 


1 + W 4) ) exp {E ~ 750) /4500] 


550 


(24) 


and 


°d = 17 


1 + exp [( E t h - E)j 12] 


- 1 


x 1 


0.21(£/145 — 1) 

1 + exp [(145 — E)j 6] 


exp (—.E/3000) (25) 


where E th is the threshold energy for the breakup reaction listed in table 2 (ref 31), and E 'is the 
kinetic energy in units of McV/amu. The low-energy behavior of equation (25) resembles that of 
the behavior in reference 28. The pickup cross section <x pickup is parameterized as 


^pickup 


48 


exp |-(E - E t h) 


1 . 7 / 3000 ] 


(26) 


and contributes to the inclusive 3 He and 2 H production cross sections. At low energies, the 
resonance 5 Li (which is not considered here) occurs. Figures 6 to 8 compare parametenzations wit 
experimental data. All cross sections are set constant above 3 GeV/amu. The energy variations 
near the thresholds and the pion production region are accurately reproduced. For the H target, 
the absorption cross section cr a b s below 80 MeV is assumed as 


CT abs = ^He + ct 3 h + a 2 H + ^pickup 


(27) 


Above 80 MeV (and below 80 MeV for A r >1), we use the energy-dependent parameterization 
of Townsend and Wilson (ref. 32). 


a abs = IOt r(3(E) [fi 4 „ e + Ra t ~ E26rf(E)] 2 


(28) 


where 


0(E) = 1 + ^ 

_ o.2 — - — I — 0.292 exp(— E/792) cos(o.228.E 0 ' 453 ) (29) 

' ' Ap Ap v 

with a normalization correction of 0.95 used for 1 H. The argument of the cosine function is in radians. 
In equation 28, the nuclear matter radii R are used. Figure 9 shows the absorption cross section ioi 
a- 1 !!. An excellent reproduction of the experimental data is seen. 

The proton and the neutron productions are expected to rise dramatically above pion produc- 
tion thresholds because two-body collisions will be predominantly inelastic, thus leading to pion 


6 




absorption in the A = 3 or A = 2 clusters. Pion production has not been treated in our theoret- 
ical considerations. However, the expected rise in proton and neutron cross sections occurs if we 
simply balance the absorption cross section with channels that do not lead to proton and neutron 
production, respectively. Figure 10 presents these results. 

The experimental data base for composite targets is extremely small. Figure 11 shows a 
comparison of equation (28) with data for the 4 He + 12 C absorption cross sections. Agreement 
is excellent, and a previous analysis (ref. 32) suggests that a similar agreement exists for other 
targets. The stripping reactions become more complicated for A r > 1 because (1) there will be 
several channels and (2) stripping or pickup to excited states of the target is a contributing factor. 
We follow the technique used by Serber in reference 33 by assuming there is a surface reaction for 

nucleon stripping on 4 He and scaling equation (26) by a\' \ Because 4 He is its own mirror nuclei, 
we ignore the coulomb effects and use equation (26) for both 3 He and 3 H production in the stripping 
reactions. Note that 3 H is not produced in the stripping reactions with *H. A slight overestimate may 
occur because we expect a small contribution for the 3 He exchange in the reaction a + l H — ► 3 He 
+ d . We also ignore any d production in the stripping process for A T > 1. Figure 12 shows results 
for 3 H and 3 He production on 12 C. The fragmentation cross section is scaled as .4^ 31 for these 
fragments. The measurements are from Webber (ref. 27), and the datum at 3.6 GcV/amu is from 
reference 34. 

Table 3 compares the parametric fits with the secondary yields for charge fragments at 
3.6 GeV/amu for several targets as measured (ref. 34). The experiment of reference 34 measured 
only peripheral events with detection angles <5°, and we expect the measurements of 1 H secondaries 
to be underestimates. The multiplicity for nucleon production from 4 He at high energies is between 
1 and 1.2 compared with a value of 2 that is assumed in existing cosmic-ray codes. For 2 H, we have 
used a scaling of Ay' 4 from our parameterization in equation (25). 

The multiplicities for proton and neutron production from deutcron projectiles are assumed as 
unity for transport calculations that ignore pickup reactions. For 3 He and 3 H projectiles, we use a 
deuteron multiplicity of 0.35 based on calculations (ref. 20), and the nucleon multiplicities are found 
by balancing the inelastic cross section. 


Secondary Energy Spectrum 

The energy spectrum of projectile fragments is determined by the internal momentum distribution 
of the projectile clusters and a weak dependence on the target caused by dynamical effects. At low 
energies, kinematical restrictions also limit the energy losses that occur. The energy spectrum of 
secondaries from light-ion fragmentation is parameterized as 

= o f N f exp [- (E - Eq + C F ) 2 /2W£] (30) 

where the width Wp and the downshift Cp are in units of MeV/amu and Eq is the beam kinetic 
energy in units of MeV/amu. In equation (30), the normalization is given by 


Np 


2 1 


n W F $ [(E -E 0 + C F )/V2W F ] - <J> [(-£(, + C F ) / (V2W f )} 
where 4>(x) is the error function 


(31) 




du 


(32) 
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and E m is the maximum kinetic energy of the fragment (if we assume a three- body final state). This 
kinetic energy is determined by the corresponding maximum momentum: 

1 /2 

p M = — ^ ({5 - [Mp + (A/p - A/p + e D ) + Mx] 2 } {5 - [(A/p - A/p + e B ) + M x - A/p] 2 }) 

VS (33) 

where y/S is the total energy in the three-body center-of-mass frame, e B is the binding energy of 
the projectile clusters, and M is the mass. The cumulative spectrum of equation (12) that results 
from equation (30) is 


(<£> { [e(r) - Eo(r') + C F ] / (72VT F ) 

}-*{ 

'[-E 0 (r') + Cp}/(V2Wp)}) 


{[E M -E 0 (r') + Cp}/(V2W F ) 

}-*{ 

{-E 0 (r')+Cp}/(V2W F )}) 


(34) 


The energy-dependent parameters Wp and Cp are described in reference 20. 

The elastically scattered energy spectrum of the alpha particle is parameterized using the Born 
term of the optical model expansion that is normalized to the coherent model results (ref. 7). A 
similar approach is followed to parameterize the quasi-elastic energy spectrum. Assuming a Gaussian 
density matrix for the target, we find from equation (22) 


da 

dE 


2m x Z oq pexp [-2 m N Z (lj - e Cl )] 8(uj - t Bx ) 


(35) 


where 6 is the unit step function and 

„ „2 (R 2 J 2 ) + B 

^(R 2 J 2) + B + R 2 t 


(36) 


The terms R n and Rj denote the matter radii of the alpha particle and target, respectively, and 
B is the slope parameter. An energy-dependent parameterization of B was given in reference 6. 
Equation (35) is expected to underestimate the spectrum at large lj because of multiple scattering 
and perhaps pion production. The cumulative spectrum that results from equation 35 is 

Fa ( r , r') = aQp {exp [—2 m x Z ( Eq — e — e#)] — exp [—2 m x Z (Eq — c^)]} (37) 

For nucleon- nucleus collisions, the energy spectrum of light target fragments is parameterized (refs. 6 
and 7) as 

-TF = , 1 i;2 ^Eexp(-£/ 2£ 0j) < 38 ) 

( 2 * 4 ) 

where E 0 j is the average energy of target fragment j. Values for E {)j are derived from the Bertini 
cascade code (ref. 35) in reference 6. For heavier projectiles, the target fragmentation cross section 
is scaled as A p 4 with the average energy kept at the nucleon value. These assumptions are expected 
to be only approximately true; however, the dominance of protons and neutrons in the GCR flux 
may diminish the importance of the accuracy of these assumptions. The most important case for 
future study will be the target fragment spectra from the incident alpha particles. The cumulative 
spectrum from equation (38) is 


Fij (r, r) = 


2 

7* 




E 


2 ’ 2 Eqj 


(39) 


where 7 (a, x) is the incomplete gamma function 

'y(a,x)= f exp (-t)t a ~ l dt 

Jo 


( 40 ) 
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The production spectra of protons and neutrons are estimated using the same mass scaling factor 
with the evaporation and cascade spectrum described in reference 6. 

Results for Light-Ion Tr ansport 

Calculations are now discussed for light-ion production using a solar minimum spectrum from 
1977 (ref. 36). Solar modulation occurs over a 22-year cycle with successive solar minimum at 
approximate 11-year intervals in which maximum particle intensities occur. The 1977 spectrum 
represents a maximum in this alternating sequence. An estimate of a primary 3 He component is also 
made. 

Figure 13 shows the energy spectra for a, h. t , and d particles at various shielding depths in 
a liquid hydrogen shield. These calculations neglect any primary d or h component. The d and 
h spectra show two local maxima: the high-energy one at ~200 McV/amu (which corresponds to 
the peak in the primary a spectrum,) and the low-energy one below 100 MeV/amu (which is due 
to the stripping reaction a -b p — ► 3 He + d ). The fragmentation spectra that are shown are 
narrow in relation to the energy grid that is normally used in the BRYNTRN code (ref. 7). The 
narrowness of the spectra leads to an inefficient numerical procedure for particles above several 
GeV/amu (illustrated in fig. 14) in which the dcuteron spectrum at 10 g/cm 2 is plotted for several 
grid-point values. Integrated quantities such as dose and dose equivalent converge quickly, and 
only approximately 90 points are needed for an accuracy of a few percent (as shown in table 4 for 
liquid hydrogen shields) and about 60 points for aluminum shields. Note that the target fragments 
contribute for the aluminum shields. 

Figure 15 gives the calculations for water shields, and figure 16 shows the calculations for 
aluminum shields. The added source of light ions from target fragmentation is now seen at low 
energies. For aluminum shields, the target fragments dominate over the projectile fragments d, t, and 
h in the total number of particles and in the expected biological effect, because the target fragments 
have lower velocities. The distribution between projectile and target effects will be much closer 
when a primary H and 3 He component is included. Approximately 10 percent of the primary He 
spectrum is the 3 He isotope with the 2 H isotopes similar in abundance to 3 He (ref. 37). Table 5 shows 
the contributions in dose and dose equivalent with a breakdown in source, either from projectiles 
or projectile fragments (indicated as PF) or from target fragmentation (indicated as TF), given for 
aluminum shields. 

Figures 17(a) and 17(b) give an estimate of the importance of a primary 3 He abundance spectrum 
in transport through liquid hydrogen and aluminum shields. We assume that the 3 He is 10 percent 
of the total helium abundance below 100 MeV/amu and that it smoothly rises to 15 percent at high 
energies (ref. 37). A comparison of figures 13 and 16 with figures 17(a) and 17(b), respectively, 
shows approximately a factor of 2 increase in 3 He at moderate shields; this increase is due to 
the primary component. Table 6 gives a comparison of the dose when 3 He is considered with and 
without the primary component for the 3 He isotopes. Including the primary 3 He leads to comparable 
contributions for the projectile- like and the target fragment 3 He dose equivalent. 

Figures 18 and 19 present a breakdown in the neutron flux from primary protons, alpha particles, 
and a combination of all heavier elements versus the shielding depth for water and aluminum shields. 
These figures demonstrate that alpha particles produce an important contribution to the total 
GCR neutron production and that their relative importance increases for lower mass shields. The 
contribution for neutrons produced from target fragmentation for alpha-induced reactions carries a 
large uncertainty because of the simple projectile mass scaling to the proton-induced spectra that we 
have used in our calculations. This uncertainty suggests a need to improve modeling for secondary 
neutron and proton and deuteron production for alpha-induced target fragmentation. Comparisons 
with experimental data (refs. 38 and 39) for monoenergetic alpha-particle beams could also be made. 

Concluding Remarks 

The space radiation transport code BRYNTRN was modified to transport cosmic-ray helium 
nuclei and their secondaries. A data base for light-ion fragmentation was discussed and implemented 
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into the transport code. Light ions produced in target fragmentation events were also transported, 
which will allow their effects on linear or lineal energy spectra to be evaluated and considered for 
further nuclear collisions by the helium particles. The contribution of neutron production by galactic 
cosmic ray (GCR) alpha particles was estimated to be significant, thus indicating that improvements 
in target fragmentation cross sections for light ions should be made in the future. The finite widths 
of the secondary spectrum from light-ion fragmentation were shown to lead to a slow convergence in 
the transport algorithms above approximately 5 GeV/amu, thus suggesting that numerical methods 
should be developed for improving this convergence. Measurements of the nucleon component in 
these events and estimates of stripping and pickup cross sections above and below pion-production 
thresholds are the most important information for validating the light-ion data base. The extension 
of the cosmic-rav space radiation transport code for isotope transport allows a more detailed study 
of the isotopic composition of the primary GCR. 


NASA Langley Research Center 

Hampton, VA 23681-0001 

July 13, 1993 
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Tabic 1. Comparisons of Calculations With Experiments for 
A = 3 Fragment Production From 4 He 

[Calculations in parentheses; experimental data from Webber (ref. 27)] 


(a) « + 12 C — > Af 


Tj a ^, MeV/amu 

<7,3 h , mb 

a 3 H( . ) m b 

203.3 

93.1 ± 9.3 (77.2) 

60.4 ± 6 (79.3) 

377.1 

79 ± 7.9 (59.9) 

66.9 ± 6.7 (60.9) 

519.9 

(62T) 1 

69.4 ± 6.9 (59.8) 


(b) a + ! H ^ A f 


Ti a ^, MeV/amu 

<73 H( .’ mb 

377.1 

26.3 ± 2.6 (19.5) 

519.9 

26.4 ± 2.6 (20.8) 

1025.0 

24.1 ± 1.9 (22.5) 


Table 2. Thresholds and Q - Values for p + 4 He 


Reaction 

<5, MeV 

(a) 

Threshold, MeV 

4 He(p, d) 3 Hc 

-18.354 

22.94 

4 He(p, 2p) 3 He 

-19.815 

24.77 

4 He(p, jm^Hc 

-20.578 

25.72 

4 He(p,/?d) 2 H 

-23.848 

29.81 

4 Hc(p, pd) 2 U 

-26.072 

32.59 

4 He(T? 7 pprm) 1 H 

-28.297 

35.37 


a Data taken from reference 31. 


Table 3. Comparison of Experimental Fragmentation Cross Sections for 
4 He Projectiles at 3.6 GeV/amu With Model Fits 


[Model fits in parentheses] 


Fragment 

Values of cr p, mb, for target s a of | 

Li 

C 

A1 

Cu 1 

'H 

166 ± 13 (536.6) 

227 ± 20 (592.0) 

319 ± 34 (823.9) 

417 ± 45 (1294.9) 

2 H 

84 ± 15 (68.2) 

91 ± 27 (91.2) 

113 ± 38 (128.2) 

159 ± 45 (184.2) 

3 h 

47 ± 5 (52.7) 

58 ± 9 (65.4) 

73 ± 20 (84.1) 

95 ± 14 (109.9) 

3 He 

48 ± 5 (48.1) 

49 ± 8 (59.6) 

70 ± 15 (76.7) 

95 ± 20 (100.2) 


Data taken from experiments in reference 27 which measured only particles in peripheral events (9 < 5°). 
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Table 5. Dose and Dose Equivalent From Light Ions in Aluminum Shields 
[Dose equivalent in parentheses] 



PF" 

TF b 

Total 

Particles 

D, cGy ( H , cSv) 

D, cGy (H, cSv) 

D, cGy (H. cSv) 

3 g/cm 2 

d 

0.011 

(0.012) 

0.342 

(1.345) 

0.353 

(1.357) 

t 

.013 

(.015) 

.082 

(.446) 

.095 

(.471) 

h 

.016 

(.026) 

.030 

(.458) 

.046 

(.484) 

a 

3.310 

(4.055) 

0.265 

(4.703) 

3.575 

(8.758) 

10 g/cm 2 

d 

0.034 

(0.035) 

0.361 

(1.465) 

0.395 

(1.500) 

t 

.038 

(0.043) 

.092 

(.508) 

.126 

(0.551) 

h 

.041 

(0.060) 

.033 

(.499) 

.074 

(0.559) 

a 

2.890 

(3.533) 

0.289 

(5.207) 

3.179 

(8.740) 

40 g/cm 2 

d 

0.144 

(0.147) 

0.342 

(1.616) 

0.486 

(1.763) 

t 

.086 

(.087) 

.097 

(.614) 

.183 

(.701) 

h 

.068 

(.079) 

.051 

(.598) 

.119 

(.677) 

a 

.071 

(.081) 

1.897 

(8.035) 

1.968 

(8.116) 


a PF represents contributions from projectile fragments. 
,J TF represents contributions caused by target fragments. 


Table 6. Dose and Dose- Equivalent Contributions From 3 He in Shielding Materials 

[Dose equivalent in parentheses] 

(a) Liquid hydrogen shields 


x, g/cm 2 

With primary 3 He 

Without primary 3 He 


D, cGy (H, 

cSv) 

D, cGy (H, cSv) 

3 

0.500 

(0.577) 

0.107 (0.141) 

5 

.485 

(.583) 

.144 (.182) 

10 

.424 

(.485) 

.183 (.219) 

20 

.286 

(.325) 

.162 (.189) 

40 

.103 

(.117) 

.068 (.079) 


(b) Aluminum shields 


x, g/cm 2 

With primary 3 He 

Without 

primary 3 He 




D, cGy ( H , cSv) 1 

D, cGy (H, cSv) 


3 

0.515 (1.026) 

0.046 

(0.484) 



5 

.511 (1.038) 

.057 

(.523) 



10 

.495 (l- 051 ) 

.074 

(.559) 



20 

.459 (1 045) 

.095 

(.662) 



40 

.390 (1.00) 

.119 

(.677) 
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-f, (GeV/c) 2 

Figure 1. Inclusive 4 He + 12 C scattering distributions at 3.6 GeV/amu versus invariant momentum transfer. 



-t, (GeV/c) 2 

Figure 2. Inclusive 4 He 4- 27 A1 scattering distributions at 3.6 GeV/amu versus invariant momentum transfer. 
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Figure 3. Angular distribution for a + 1 II — > ,5 He at 1.02 GeV/amu. 




Figure 4. Comparison of calculations of longitudinal momentum distribution for 3 H production in 4 He + 12 C 
collisions at 1.9 GeV/amu with experiment (ref. 23). 



Figure 5. Comparison of calculations of transverse momentum distribution for 3 H production in 4 He + 12 C 
collisions at 0.385 GeV/amu with experimental data. 
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Figure 6. Comparison of parametric model for ^He production in a 4- *H collisions with experimental data. 



Figure 7. Comparison of parametric model for 3 H production in a + ! H collisions with experimental data. 
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Figure 8. Comparison of parametric model 2 H production in a + 'll collisions with experimental data. 


Calculations 



Figure 9. Comparison of parametric model for ol + absorption cross section with experimental data. 


20 


Proton 

Neutron 


iT] b* MeV/amu 

Figure 10. Parametric model predictions for proton production and neutron production in a + 1 H collisions, 

800 pr Calculation 

\ t Experimental with error 

700 

-\ ' 

- 


\ l 


-D 

cS 

t> 

\ 

i 

500 

— 

^ T 

400 

1 1 1 1 1 L 

i 

1 1 1 1 1 1 1 L 1 t t 1 1 I | L 1 _i t 1 1 1 


21 




$ (experimental with error) 
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Model 



Figure 12. Comparison of fits for 3 He and 3 He production cross sections in 4 He+ 12 C collisions with experimental 
data. 
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Figure 14. Comparison of deuteron flux spectra in 10 g/cm 2 of H 2 shields for several grid values. 
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Figure 16. Calculations of light-i 
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Figure 18. Secondary neutron flux contributed versus depth in water shield. 



Figure 19. Secondary neutron flux contributed versus depth for aluminum shield. 
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